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61  Summary 

62  During  routine  screening  for  endemic  Burkholderia  pseudomallei  from  water  wells  in  northern 

63  Australia,  Gram-negative  bacteria  (strains  MSMB43  T,  MSMB121,  and  MSMB122)  with  a 

64  similar  morphology  and  biochemical  pattern  to  B.  pseudomallei  and  B.  thailandensis  were  co- 

65  isolated  with  B.  pseudomallei  on  Ashdown’s  selective  agar.  To  determine  the  exact  taxonomic 

66  position  of  these  strains  and  to  distinguish  them  from  B.  pseudomallei  and  B.  thailandensis,  they 

67  were  subjected  to  a  series  of  phenotypic  and  molecular  analyses.  Biochemical  and  fatty  acid 

68  methyl  esters  analysis  was  unable  to  distinguish  B.  humptydooensis  sp.  nov.  from  closely  related 

69  species.  In  MALDI-TOF  analysis,  all  isolates  grouped  together  in  a  cluster  separate  from  other 

70  Burkholderia  spp.  16S  rRNA  and  recA  sequence  analysis  demonstrated  phylogenetic  placement 

71  for  B.  humptydooensis  sp.  nov.  in  a  novel  clade  within  the  B.  pseudomallei  group.  MLST 

72  analysis  of  the  three  isolates  in  comparison  with  MLST  data  from  3,340  B.  pseudomallei  strains 

73  and  related  taxa  revealed  a  new  sequence  type  (ST3 18).  Genome  to  genome  distance  calculations 

74  and  average  nucleotide  identity  of  all  isolates  to  both  B.  thailandensis  and  B.  pseudomallei, 

75  based  on  whole  genome  sequences,  also  confirmed  B.  humptydooensis  sp.  nov.  as  a  novel 

76  Burkholderia  species  within  the  pseudomallei  complex.  Molecular  analyses  clearly  demonstrate 

77  that  strains  MSMB431,  MSMB121,  and  MSMB122  belong  to  a  novel  Burkholderia  species  for 

78  which  the  name  Burkholderia  humptydooensis  sp.  nov.  is  proposed  with  the  type  strain 

79  MSMB431  (=ATCCXXXXXt;  LMGXXXX,  CP013380-CP013382). 

80  Burkholderia  species  are  abundant  and  occupy  diverse  ecological  niches,  including  soil, 

81  plants,  animals,  and  humans.  Probably  the  most  diverse  and  environmentally  adaptable  plant- 

82  associated  bacteria  also  belong  to  the  genus  Burkholderia  (Compant  et  ah,  2008).  Many  species 

83  of  Burkholderia  have  been  described  since  the  discovery  of  B.  cepacia  by  W.  H.  Burkholder  in 
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84  1949  as  the  cause  of  onion  rot  (Burkholder,  1950),  and  later  recognized  as  a  human  pathogen. 

85  Currently,  there  are  more  than  90  identified  species  in  this  genus  (Vial  et  ah,  2011,  Gu  et  ah, 

86  2015).  There  has  been  a  proposal  to  divide  the  species  into  two  genera,  one  of  which  would 

87  retain  the  Burkholderia  name  and  the  other  Paraburkholderia  gen.  nov.  (Sawana  et  ah,  2014).  At 

88  least  17  closely-related  species  belong  to  the  Burkholderia  cepacia  complex  (BCC),  including 

89  nine  well-known  species  comprising  the  genomovars  I-IX,  respectively:  B.  cepacia,  B. 

90  multivorans,  B.  cenocepacia,  B.  stabilis,  B.  vietnamiensis,  B.  dolosa,  B.  ambifaria,  B.  anthina,  B. 

91  pyrrocinia  (Coenye  and  Vandamme,  2003);  seven  recently  identified  species:  B.  latens,  B. 

92  diffusa,  B.  arboris,  B.  seminalis ,  B.  metallica,  B.  contaminans,  and  B.  lata  (Vanlaere  et  ah,  2009, 

93  Vanlaere  et  al.,  2008);  and  B.  ubonensis  (Yabuuchi  et  al.,  2000).  Most  of  these  species  reside  in 

94  soils  and  are  associated  with  plants. 

95  Notably,  there  are  two  Burkholderia  species  that  can  cause  severe  human  and  animal 

96  diseases:  B.  pseudomallei  and  B.  mallei,  the  causative  agents  of  melioidosis  and  glanders, 

97  respectively.  B.  pseudomallei  is  a  major  cause  of  community  acquired  sepsis  in  Northeast 

98  Thailand  and  northern  Australia  (Cheng  and  Currie,  2005).  Due  to  the  concerns  of  their  potential 

99  use  as  weapons  of  mass  destruction,  federal  health  agencies  in  the  United  States  have  recently 

100  classified  these  species  as  Tier  1  (Top  Tier)  disease  agents  (Butler,  2012).  It  has  been  well 

101  established  that  B.  mallei  is  a  clone  of  B.  pseudomallei  that  became  a  host-adapted  pathogen  in 

102  equines  resulting  in  massive  genome  reduction  (Godoy  et  al.,  2003).  Genetically,  both  species 

103  are  members  of  the  pseudomallei  phylogenetic  group  or  complex 

104  (http://www.ncbi.nhn.nih.gov/Taxonomy/).  Three  additional  closely  related  species  have  been 

105  identified  so  far  in  this  group,  including  B.  thailandensis  (Brett  et  al.,  1997),  B.  oklahomensis 

106  (Glass  et  al.,  2006b),  and  a  newly  identified  B.  thailandensis -like  species  (Currie,  2015,  Gee  et 
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107  al.,  2008).  These  closely-related  species  are  soil  saprophytes  and  are  considered  non-pathogenic, 

108  although  a  few  strains  of  B.  thailandensis  and  B.  oklahomensis  have  been  described  as  causing 

109  clinical  infection  in  humans  (Glass  et  al.,  2006a,  Glass  et  ah,  2006b). 

110 

111  Strain  MSMB431  (proposed  type  strain)  was  isolated  from  a  water  sample  from  an 

112  automated  water  bore  (well)  collected  in  1995  and  examined  for  B.  pseudomallei  in  the  Northern 

113  Territory  (NT)  of  Australia  (Gee  et  al.,  2008).  It  was  initially  thought  to  be  B.  thailandensis 

114  because  this  strain  was  able  to  assimilate  arabinose  as  a  sole  carbon  source,  which  is  a  trait  used 

115  to  discriminate  B.  thailandensis  from  B.  pseudomallei  (arabinose  negative)  (Gee  et  al.,  2008). 

116  The  bore  from  which  MSMB43T  was  discovered  is  located  in  Humpty  Doo,  a  region  of  rural 

117  properties  outside  the  capital  of  the  NT,  Darwin.  The  Top  End  of  the  NT  has  a  high  incidence 

118  rate  of  melioidosis  (Parameswaran  et  al.,  2012).  Indeed  the  water  sample  from  which  MSMB431 

119  was  recovered  also  yielded  B.  pseudomallei.  An  additional  two  strains  (MSMB121  and 

120  MSMB 122)  were  both  isolated  in  2007,  from  a  single  separate  bore  water  sample  within  the  NT 

121  collected  approximately  950  km  south  of  the  territory  capital,  Darwin,  resulting  in  a  910  km 

122  separation  between  the  two  sample  sites  of  MSMB121/MSMB122  andMSMB43T.  To  date,  B. 

123  humptydooensis  sp.  nov.  has  not  been  identified  outside  the  NT  and  it  has  not  been  isolated  from 

124  any  clinical  specimens  from  patients  within  the  NT. 

125  The  specific  epithet  “ humptydooensis ”  given  to  this  new  species  was  adopted  from  the 

126  location  name  “Humpty  Doo,”  the  location  where  this  new  species  was  first  discovered.  The 

127  name  “humptydooensis”  has  been  used  previously  for  a  new  species  of  beetle  ( Australoxenella 

128  humptydooensis)  from  this  same  geographic  region  (Howden  and  Storey,  1992). 

129 
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130  As  described  previously  (Gee  et  al.,  2008),  MSMB431  did  not  grow  while  incubated  at  a 

131  temperature  greater  than  42°C  and  also  produced  little  or  no  gas  from  nitrate.  All  three  strains 

132  were  grown  at  8,  25,  37,  42,  and  45°C  for  24,  48,  72,  and  144  h  on  Columbia  blood  agar, 

133  MacConkey  agar,  Ashdown’s  selective  agar,  and  Standard  I  with  and  without  supplementary 

134  CO2.  On  Columbia  blood  agar  smooth  and  creamy-white  colonies  were  observed  after  24  h, 

135  whereas  red,  convex  and  small  (l-2mm)  colonies  were  observed  on  MacConkey  after  48h.  Dry 

136  and  wrinkled  colonies  were  observed  on  Ashdown’s  agar  after  72  h  of  growth  (Fig.  1),  similar  to 

137  the  appearance  of  B.  pseudomallei,  while  slimy-confluent  honey-like  growth  appeared  on 

138  Standard  I  after  48h  (Fig.  1).  Bacterial  growth  was  visible  on  all  media  after  incubation  at  25- 

139  42°C  for  at  least  24  h  with  best  growth  observed  on  Columbia  blood  agar.  No  growth  was 

140  observed  at  8°C  and  45°C.  The  optimal  temperatures  for  growth  were  between  28-37°C 

141  aerobically. 

142  Cell  morphology  was  examined  using  a  Zeiss  light  microscope  at  lOOOx  magnification 

143  with  cells  grown  for  2  days  at  37  °C.  All  strains  showed  a  Gram-negative  bipolar  staining 

144  appearing  as  rods  of  2-3  pm  in  length  and  0.4-0. 8  pm  in  diameter.  All  strains  were  motile  in 

145  semi-solid  media. 

146  Biochemical  data  were  obtained  for  all  three  strains  of  B.  humptydooensis  sp.  nov. 

147  (MSMB431,  MSMB 12 1 ,  and  MSMB 122)  and  compared  to  the  type  strains  of  B.  pseudomallei 

148  (K96243)  (Wuthiekanun  et  al.,  1996)  and  B.  thailandensis  (E264)  using  the  API  NE  and  API 

149  Zym  (bioMerieux)  systems  according  to  the  manufacturer’s  instructions.  Differentiation  of  B. 

150  humptydooensis  sp.  nov  from  B.  thailandensis  was  not  possible  based  on  the  biochemical 

151  reactions  (Table  1). 
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152  Matrix-assisted  laser-desorption/ionization  time-of-flight  mass  spectrometry  (MALDI- 

153  TOF  MS)  was  performed  on  all  three  B.  humptydooensis  sp.  nov  strains  and  they  clustered  with 

154  other  members  of  the  B.  pseudomallei  complex  (see  SI  doc  and  Fig.  SI).  Fatty  acid  methyl  esters 

155  analysis  was  unable  to  distinguish  among  the  fatty  acid  profiles  from  the  three  B. 

156  humptydooensis  sp.  nov.  strains  and  its  closely  related  species  (five  B.  thailandensis,  two  B. 

157  oklahomensis,  and  three  B.  ubonensis  strains)  (see  SI  doc  and  Fig.  S2). 

158  Minimum  inhibitory  concentrations  (MICs)  were  determined  by  broth  microdilution 

159  method  using  commercially  available  CE-certified  MICRONAUT-S  96-well  microtitre  plates 

160  (Merlin,  Bomheim-Hersel,  Gennany)  containing  a  2-fold  serial  dilution  of  the  following 

161  antibiotics:  amoxicillin/clavulanic  acid  (0.5-64/0.25-32  mg/L),  ceftazidime  (0.5-64  mg/L), 

162  imipenem  (0.25-32  mg/L),  rifampicin  (0.0625-8  mg/L),  chloramphenicol  (0.5-64  mg/L), 

163  trimethoprim/sulfamethoxazole  (0.25-32/4.75-608  mg/L),  streptomycin  (0.25-32  mg/L), 

164  gentamicin  (0.25-32  mg/L),  doxycycline  (0.25-32  mg/L),  tigecycline  (0.03 125-4  mg/L), 

165  ciprofloxacin  (0.03 125-4  mg/L),  and  levofloxacin  (0.0625-4  mg/L).  One  well  without  antibiotic 

166  was  used  as  a  growth  control.  All  plates  containing  the  lyophilized  antimicrobial  substances  were 

167  stored  at  room  temperature  until  use. 

168  Testing  conditions  were  in  accordance  with  the  current  Clinical  &  Laboratory  Standards 

169  Institute  (CLSI)  recommendations  for  B.  pseudomallei  (CLSI  M45-A2  Volume  30  No.  18). 

170  Single  colonies  of  MSMB431,  MSMB121,  and  MSMB122  were  picked  from  agar  plates  and 

171  inoculated  in  physiological  saline  until  the  turbidity  matched  that  of  a  0.5  McFarland  standard. 

172  The  suspension  was  diluted  22 1 -fold  in  cation-adjusted  Mueller  Hinton  II  broth  (Becton 

173  Dickinson,  catalog  #297701).  After  incubation  for  24  h  at  37°C  in  a  5%  CO2  atmosphere 

174  bacterial  growth  was  verified  photometrically  at  a  wavelength  of  620  nm  using  a  commercial 
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175  photometer  (MERLIN,  Bomheim-Hersel,  Germany)  and  each  strain  was  tested  in  triplicate. 

176  Additionally,  a  gradient  strip  method  (Etest,  bioMerieux)  was  applied  to  investigate  a  broader 

177  range  of  antibiotic  concentrations. 

178  Adopting  the  CLSI  breakpoints  of  B.  pseudomallei,  all  strains  were  determined  to  be  in 

179  vitro  susceptible  to  ceftazidime,  imipenem,  trimethoprim/sulfamethoxazole,  and  doxycycline, 

180  whereas  a  resistance  to  amoxicillin/clavulanic  acid  was  observed.  (Table  2).  The  antimicrobial 

181  susceptibility  pattern  of  B.  humptydooensis  sp.  nov.  generally  resembles  that  of  B.  pseudomallei 

182  (Ahmad  et  al.,  2013,  Thibault  et  ah,  2004,  Trunck  et  ah,  2009,  Crowe  et  ah,  2014),  except  for  the 

183  elevated  MICs  of  amoxicillin/clavulanic  acid  (Table  2).  No  significant  differences  were  observed 

184  either  among  the  three  strains  or  between  the  two  different  susceptibility  testing  methods.  As  the 

185  maximum  concentration  of  aminoglycosides  in  the  microtiter  plates  was  32  mg/L,  high-level 

186  streptomycin  resistance  but  low-level  gentamicin  resistance  was  confirmed  using  the  Etest 

187  method  (data  not  shown). 

188  The  pathogenic  potential  of  B.  humptydooensis  sp.  nov.,  MSMB43  ,  was  investigated  in  a 

189  BALB/c  mouse  model  and  compared  to  B.  thailandensis  (type  strain  E264).  Live  culture  was 

190  grown  to  logarithmic  phase  (OD600  ~  1-0)  in  Luria-Bertani  (LB)  broth  as  described  in  (Morici  et 

191  ah,  2010).  Sterile  lxPBS  was  used  to  wash  cells  twice  before  making  dilutions  for  injecting 

192  mice.  Viability  counts  of  the  final  inocula  were  made  on  LB  agar  plates.  Six  to  eight  week  old 

193  female  BALB/c  mice  in  treatment  groups  of  5  mice  per  cage  were  used.  Lood  and  water  were 

194  provided  ad  libitum.  All  mice  in  a  single  cage  received  the  same  infectious  dose  ( B . 

195  humptydooensis  sp.  nov.:  1.05  x  1 04,  105  or  106  CPU;  B.  thailandensis :  3.4x  104,  105  or  106 

196  CPU)  via  a  single  subcutaneous  (s.c.)  injection  in  the  scruff  of  the  neck.  Mice  were  monitored 

197  daily  for  health  status.  All  mice  were  euthanized  on  day  2 1  post  injection.  This  work  was 
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198  conducted  under  approved  protocols  from  the  NAU  IACUC  (Protocol  14-0 11)  and  DOD 

199  ACURO  (HDTRAl-12-C-0066_Wagner).  Neither  species  caused  mortality  in  any  mice  when 

200  delivered  via  the  s.c.  route,  nor  did  any  mice  show  outward  signs  of  illness.  In  comparison,  s.c. 

201  infections  of  fully  virulent  B.  pseudomallei  results  in  50%  mortality  within  10  days  at  a  dose  of 

202  10  CFU  (Bames  and  Ketheesan,  2005).  It  remains  unknown  if  the  inhalation  route  increases  the 

203  pathogenicity  of  species  tested  in  the  same  way  as  B.  thailandensis  E264,  which  can  cause  high 

204  mortality  in  mice  at  doses  of  104  -  106  CFU  when  delivered  as  an  aerosol  (Morici  et  ah,  2010, 

205  West  et  ah,  2008,  Wiersinga  et  ah,  2008). 

206  16S  rRNA  and  recA  gene  sequencing  analyses  were  performed  on  three  B. 

207  humptydooensis  sp.  nov.  strains:  MSMB431,  MSMB121,  and  MSMB122  as  previously  described 

208  (Brett  et  ah,  1997,  Payne  et  ah,  2005).  From  whole  genome  analysis  of  strain  MSMB43Tusing 

209  SSU-ALIGN  (Nawrocki  et  ah,  2009)  we  detennined  that  four  rRNA  operons  are  present  on  its 

210  chromosomes  of  which  two  unique  versions  were  found  (AQ610  12930,  AQ610  01425; 

211  AQ610_21350,  AQ610_02540).  These  two  copies  of  the  16S  rRNA  genes  are  different,  which 

212  leads  to  ambiguities  in  conventional  sequencing  (see  Fig.  S3).  The  16S  rRNA  gene  sequence 

213  similarities  of  B.  humptydooensis  sp.  nov.  to  other  members  of  the  pseudomallei  complex  (B. 

214  thailandensis ,  B.  mallei,  B.  oklahomensis )  were  99%.  Phylogenetic  reconstruction  of  16S  rRNA 

215  and  recA  sequences  using  MEGA  version  6  (Tamura  et  ah,  2013)  confirmed  genetic  proximity  to 

216  the  pseudomallei  complex  but  also  determined  that  all  B.  humptydooensis  sp.  nov.  strains  formed 

217  their  own  group  within  this  complex  (see  Fig.  S3-S4). 

218  MLST  was  performed  on  all  three  B.  humptydooensis  sp.  nov.  strains  as  previously 

219  described  (Godoy  et  ah,  2003).  As  of  October  6,  2015,  a  total  of  1,429  sequence  types  (STs)  had 


220  been  identified  in  B.  pseudomallei  and  closely  related  species  by  MLST  (http://www.MLST.net). 
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221  All  three  B.  humptydooensis  sp.  nov.  strains,  MSMB431,  MSMB121,  and  MSMB122  are  ST318 

222  and  there  are  no  other  representatives  of  this  ST.  Overall,  phylogenetic  analysis  using  MLST 

223  data  supports  the  separation  of  B.  humptydooensis  sp.  nov.  from  the  other  pseudomallei  complex 

224  members,  as  described  previously  (Gee  et  ah,  2008). 

225  Two  genomes  were  sequenced  on  the  PacBio  platfonn,  resulting  in  one  finished  assembly 

226  (MSMB431)  and  one  mostly  finished  assembly  with  4  contigs  (MSMB122)  (Table  3).  Two  other 

227  genomes  that  group  with  B.  humptydooensis  sp.  nov.  are  present  in  GenBank,  see  Figure  2  ( B . 

228  thailandensis  MSMB121and  Burkholderia  spp.  MSMB43  (2002721687)  with  GenBank 

229  assembly  accession  numbers  of  GCA_000385525  and  GCA_000959325  respectively).  The 

230  MSMB431  genome  had  one  circular  contig  ~305Kb  long  that  appears  to  be  a  plasmid;  this  same 

231  sequence  is  also  present  in  the  completed  genome  of  Burkholderia  spp.  MSMB43  (2002721687). 

232  A  comparative  genomics  approach  using  LS-BSR  (Sahl  et  ah,  2014)  demonstrated  that  a  large 

233  stretch  of  the  B.  pseudomallei  K96243  genome  (BPSS1 165  -  BPSS1 184)  on  chromosome  2  is 

234  highly  conserved  (>98%  identity)  in  the  plasmid  sequence,  suggesting  a  shared  origin  for  these 

235  regions. 

236  For  the  core  genome  phylogeny,  genomes  were  aligned  against  B.  pseudomallei  K96243  with 

237  NUCmer  (Delcher  et  al.,  2003).  The  reference  genome  was  also  aligned  against  itself  to  identify 

238  duplicated  regions,  which  were  masked  from  subsequent  analyses;  the  NASP  pipeline  was  used 

239  to  wrap  these  methods  (http://tgennorth.github.io/NASP/).  A  phylogeny  was  inferred  with 

240  RAxML  v8  (Stamatakis,  2014)  on  a  large  set  (n=33 1,000)  of  concatenated  SNPs  using  a  time 

241  reversible  model  incorporating  the  Lewis  ascertainment  bias  correction.  The  phylogeny 

242  demonstrates  the  position  of  B.  humptydooensis  sp.  nov.  in  relation  to  other  clades  in  the  B. 

243  pseudomallei  complex  and  confirms  the  results  from  other  methods  (Fig.  2). 
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DNA-DNA  hybridization  (DDH)  is  the  current  gold-standard  for  bacterial  species 
delineation.  DDH  is  necessary  for  the  description  of  a  new  species  within  a  taxon  when  strains 
share  more  than  97%  of  16S  rRNA  gene  sequence  similarity  (Goebel  and  Stackebrandt,  1994, 
Tindall  et  ah,  2010).  If  DNA-DNA  relatedness  between  two  strains  is  less  than  70%,  both  strains 
are  considered  as  different  species.  However,  DDH  is  laborious  and  difficult  to  standardize  and 
inter-laboratory  reproducibility  is  relatively  low.  DDH  was  previously  performed  on  MSMB43T 
and  has  a  relative  binding  ratio  of  91%  with  B.  thailandensis  (ATCC  700388T)  with  a 
divergence  of  4%  (Gee  et  ah,  2008).  Because  of  the  drawbacks  of  conventional  DDH  and  the 
rapid  progress  in  genome  sequencing  techniques,  various  in  silico  algorithms  for  calculating 
genome-to-genome  similarities  or  distances  have  been  developed.  In  addition  to  the  commonly 
used  average  nucleotide  identity  (ANI)  (Konstantinidis  and  Tiedje,  2005),  recently  a  highly 
reliable  estimator  for  the  relatedness  of  genomes  was  developed  by  Jan  P.  Meier-Kolthoff  and 
colleagues  (Meier-Kolthoff  et  ah,  2013).  Genome  to  genome  distance  calculation  (GGDC) 
produces  digital  DDH  values  that  correlate  well  with  values  obtained  by  conventional  DDH, 
which  is  of  utmost  importance  for  compatibility  with  the  current  species  concept,  and  also 
provides  confidence-interval  estimation.  Values  reported  in  Table  4  were  generated  using 
formula  2,  which  summed  the  identities  found  in  high-scoring  segment  pairs  (HSP)  and  then 
were  divided  by  the  overall  HSP  length.  The  GGDC  service  is  available  from  the  German 
Collection  of  Microorganisms  and  Cell  Cultures  homepage  (http:// ggdc . dsmz.de/ distcalc2  .php). 
PacBio  assemblies  were  used  to  determine  the  distances  among  B.  humptydooensis  sp.  nov.  and 
other  closely  related  species,  B.  humptydooensis  sp.  nov.  strains  MSMB431,  MSMB121,  and 
MSMB122.  The  genomes  were  subjected  to  GGDC  analysis  and  compared  to  the  available 
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267  genome  sequences  of  B.  pseudomallei,  B.  oklahomensis  and  B.  thailandensis  reference  strains 

268  (K96243,  C6786t,  and  E264,  respectively).  For  comparison,  the  average  nucleotide  identity 

269  (ANI)  values  were  also  calculated  for  all  reference  sequences  with  JSpecies  (Richter  and 

270  Rossello-Mora,  2009);  the  authors  of  JSpecies  detennined  that  ANI  values  <95%  indicate 

271  separate  species. 

272  Among  the  three  tested  B.  humptydooensis  sp.  nov.  genomes,  the  calculated  GGDC  and 

273  ANI  values  were  in  the  range  of  93-99%  and  98-99%,  respectively  (Table  4),  indicating  that  all 

274  of  these  tested  strains  belong  to  a  single  species,  including  the  proposed  B.  humptydooensis  sp. 

275  nov.  type  strain  MSMB431.  As  expected  from  WGS  analyses,  strain  MSMB431  had  a  slightly 

276  lower  GGDC  similarity  (93%)  to  the  other  two  B.  humptydooensis  strains,  which  were 
Til  approximately  97%  similar.  Further  GGDC  analysis  (Table  4)  has  determined  that  the 

278  similarities  of  all  B.  humptydooensis  strains  compared  to  all  other  tested  Burkholderia  species  in 

279  the  B.  pseudomallei  complex  was  less  than  70%,  with  the  highest  similarity  detected  to  B. 

280  thailandensis  (51.1%  (±  3.2%)).  This  confirms  that  the  three  tested  strains  are  not  B. 

281  thailandensis ,  but  rather  a  distinct  species.  The  GGDC  similarity  between  B.  mallei  and  B. 

282  pseudomallei  was  92.5%  which  confirms  previous  conventional  DDH  results  and  demonstrates 

283  that,  from  a  strict  taxonomic  point  of  view,  they  belong  to  a  single  species. 

284  Description  of  Burkholderia  humptydooensis  sp.  nov. 

285  Burkholderia  humptydooensis  sp.  nov.  (hump.ty.doo.en’sis.  F.  gen.  adj.  humptydooensis, 

286  pertaining  to  Humpty  Doo,  a  small  town  in  Northern  Territory  of  Australia,  where  the  first 

287  member  of  this  species  was  isolated). 

288  Bacilli,  0.4-0. 8  pm  in  diameter  and  2-3  pm  in  length,  arranged  individually  or  in  irregular 

289  clusters.  The  organism  is  Gram-negative  with  bipolar  staining,  motile,  and  non-spore  forming. 
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290  Growth  is  observed  in  a  temperature  range  of  25-42°C  within  24-48  h  on  various  standard  solid 

291  media.  Within  24  h,  small  colonies  (0.5-1  mm)  are  formed  on  non-selective  media  (Columbia 

292  blood  and  Standard-I)  and  after  48  h  also  on  selective  media  (Ashdown’s,  MaConkey).  Best 

293  growth  occurs  at  28-37°C  after  >  24  h.  Colonies  become  confluent  and  honey-like  in  appearance 

294  on  glycerol-  containing  medium  (ST-I)  after  48  h.  On  Ashdown’s  selective  agar,  highly  wrinkled 

295  purple  colonies  are  observed  within  >  48  h,  resembling  the  growth  of  B.  pseudomallei. 

296  Positive  assimilation  (API  20NE)  was  found  for  D-glucose,  L-arabinose,  D-mannose,  D- 

297  mannitol,  N-acetyl-glucosamine,  potassium  gluconate,  capric  acid,  adipic  acid,  malic  acid, 

298  trisodium  citrate,  and  phenylacetic  acid  whilenegative  for  D-maltose.  Aesculin  and  gelatin  are 

299  hydrolyzed.  Variable  reactions  are  L-arginine  and  4-Nitrophenyl-P  D-galactopyranoside  (PNPG). 

300  Positive  (API  ZYM)  for  alkaline  phosphatase,  esterase,  esterase  lipase,  lipase,  leucine 

301  arylamidase,  acidic  phosphatase,  and  naphthol-AS-BI-phosphohydrolase.  Enzymes  absent  on 

302  API  ZYM  are  valine  arylamidase,  cystin  arylamidase,  trypsin,  a-chymotrypsin,  a-  and  [!- 

303  galactosidase,  p-glucuronidase,  a-  and  p-glucosidase,  N-acetyl-P-glucosaminidase,  a- 

304  mannosidase,  and  a-frucosidase.  This  species  is  aerobic,  catalase  and  oxidase-positive,  urease 

305  and  indole  negative.  Nitrate  and  nitrite  are  reduced  (with  no  gas  fonnation  from  nitrite)  and  no 

306  production  of  H2S. 

307  B.  humptydooensis  sp.  nov.  strains  are  resistant  to  aminoglycosides  and 

308  amoxicillin/clavulanic  acid,  but  susceptible  to  trimethoprim/sulfmethoxazole,  doxycycline, 

309  imipenem,  and  ceftazidime.  All  B.  humptydooensis  sp.  nov.  strains  are  seroreactive  with  sera 

310  from  melioidosis  patients  who  were  infected  with  B.  pseudomallei  serotype  B  strains.  All  strains 

311  produced  O-antigen  ladder  type  B2,  except  that  strain  MSMB431  produced  a  novel  O-antigen 

312  ladder  type  (Stone  et  ah,  2012).  The  type  strain,  MSMB431,  has  been  previously  referred  to  as  B. 
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313  thailandensis- like  species  in  multiple  studies  (Currie,  2015,  Gee  et  al.,  2008).  MSMB43 1  was 

314  isolated  in  1995  from  an  automated  water  well  (bore)  in  Humpty  Doo,  Australia.  Of  interest,  B. 

315  pseudomallei  was  also  cultured  from  that  same  water  sample.  B.  humpty  doo  ensis  sp.  nov.,  like  B. 

316  thailandensis,  is  non-virulent  in  mice.  In  addition,  MSMB431  is  known  to  produce 

317  Thailanstatins,  which  possess  antiproliferative  activities  in  representative  human  cancer  cell  lines 

318  (Liu  et  al.,  2013).  The  type  strain  MSMB431  has  been  deposited  in  the  American  Type  Culture 

319  Collection  as  ATCC  XXXXX  and  the  Belgian  Co-ordinated  Collections  of  Micro-organisms  as 

320  LMGXXXX. 
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328  Figure  1.  Colony  morphology  of  B.  humptvdooensis  sp.  nov.  MSMB431 

329  (a)  Culture  grown  on  Ashdown’s  agar  and  on  (b)  Standard  I. 


333 

334 

335 

336 


Page  15  of  23 


337  Figure  2.  Core  genome  phylogeny  of  B.  humptydooensis  sp.  nov. 

338  SNPs  from  the  comparison  of  four  B.  humptydooensis  sp.  nov.  genomes  and  representatives  of 

339  the  other  closely  related  species  were  used  to  reconstruct  the  phylogenetic  relationships. 

340  Genomes  from  this  study  in  bold  and  assembly  numbers  in  parentheses.  Numbers  at  nodes 

341  indicate  bootstrap  support  values.  Collapsed  nodes  shown  in  gray. 


342 

343 

344 
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345  Table  1.  Phenotypic  characteristics  of  B.  humptydooensis  sp.  nov.  and  closely  related  species 

346  within  the  B.  pseudomallei  group. 

347  Strains:  Bp,  Burkholderia  pseudomallei  (K96243)*;  Bt,  B.  thailandensis  (E264);  Bh-43,  B. 

348  humptydooensis  sp.  nov.  (MSMB431);  Bh-121.  B.  humptydooensis  sp.  nov.  (MSMB121);  Bh- 

349  122,  B.  humptydooensis  sp.  nov.  (MSMB122).  +,  Positive  reaction;  -,  negative  reaction.  All 

350  strains  were  positive  for  nitrate,  gelatin,  glucose,  mannose,  mannitol,  N-acetylglucosamine, 

351  gluconate,  caprate,  adipate,  malate,  citrate,  and  phenylacetate.  All  strains  were  negative  for 

352  glucose  (acidification)  and  urea. 


Characteristic 

^  * 

Bp 

Bt 

Bh-43 

Bh-122 

Bh-121 

Tryptophan 

+ 

- 

- 

- 

- 

Arginine 

+ 

- 

+ 

- 

- 

Aesculin 

- 

+ 

+ 

+ 

+ 

PNPG 

- 

- 

- 

- 

+ 

Assimilation  of: 

Arabinose 

- 

+ 

+ 

+ 

+ 

Maltose 

- 

+ 

- 

- 

- 

353  *  Data  obtained  from  Wuthiekanun  et  ah,  (1996). 

354 
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355  Table  2.  Summary  of  MIC’s  determined  in  triplicate  by  broth  microdilution  method. 

356  Concentrations  are  given  in  mg/L. 


Antimicrobial  substance 

B. 

pseudomallei 

MIC,0 

MSMB43t 

MSMB121 

MSMB122 

Amoxicillin/clavulanic  acid 

8° 

32/16* 

32/16* 

32/16* 

Ceftazidime 

4a 

4 

2 

4 

tmipenem 

r 

0.5 

0.5 

0.5 

Rifampicin 

32° 

>8 

>8 

8 

Chloramphenicol 

16° 

8 

4 

4 

Trimethoprim/sulfamethoxazole 

1.5e 

0.5/9. 5 

<0.25/4.75 

<0.25/4.75 

Streptomycin 

n/a 

>32 

>32 

>32 

Gentamicin 

64° 

32 

32 

>32 

Doxycycline 

la 

1 

1 

1 

Tigecycline 

4h 

2 

4 

4 

Ciprofloxacin 

4“ 

0.5 

1 

1 

Levofloxacin 

4a 

0.5 

1 

1 

357 

sjs 

358  Resistance  was  observed  based  upon  the  CLSI  breakpoints  of  B.  pseudomallei. 

359  Data  obtained  from  a  Thibault  et  al.,  (2004)  b  Ahmad  et  al.,  (2013)  c  Crowe  et  al.,  (2014). 

360 
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361  Table  3.  Whole  genome  data  for  B.  pseudomallei  group. 


Strain/species 

GC  content 

Genome  size 

#  CDS 

MSMB431 

B.  humptydooensis 

67.1% 

7.3*  Mb 

6324 

MSMB121 

B.  humptydooensis 

67.5% 

6.7  Mb 

5795 

MSMB122 

B.  humptydooensis 

67.5% 

6.8  Mb 

5845 

E264 

B.  thailandensis 

67.6% 

6.7  Mb 

5652 

C6786 

B.  oklahomensis 

66.9% 

7.1  Mb 

6097 

K96243 

B.  pseudomallei 

68.1% 

7.2  Mb 

5948 

ATCC23344 

B.  mallei 

68.5% 

5.8  Mb 

5006 

Two  chromosomes  present  for  all  genomes. 
*  One  plasmid  present. 

CDS  =  coding  DNA  sequences 


362 

363 

364 

365 

366 

367 
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368  Table  4.  Genome  to  genome  distance  calculation  (GGDC)  in  the  bottom  half  matrix  (with 

369  confidence  intervals)  and  average  nucleotide  identity  (ANI)  in  the  top  half  matrix  that 

370  demonstrate  whole  genome  sequence  similarities  of  B.  humptydooensis  sp.  nov.  strains  and  other 

371  species  in  the  pseudomallei  complex.  Values  highlighted  yellow  represent  values  above  the 

372  similarity  threshold  that  defines  members  of  the  same  species. 

373 


MSMB43 

B.  humptydooensis 

MSMB121 

B.  humptydooensis 

MSMB122 

B.  humptydooensis 

E264 

B.  thailandensis 

C6786 

B.  oklahomensis 

K96243 

B.  pseudomallei 

ATCC23344 

B.  mallei 

MSMB43 

B.  humptydooensis 

100 

98.81 

98.98 

93.56 

91.80 

91.17 

93.22 

MSMB121 

B.  humptydooensis 

93.1  ±2.3 

100 

99.84 

93.59 

91.80 

93.39 

93.26 

MS  MB  122 

B.  humptydooensis 

93.1±2.3 

97.310.02 

100 

93.60 

91.86 

93.23 

93.25 

E264 

B.  thailandensis 

51.113.2 

51.413.2 

51.313.2 

100 

91.48 

93.01 

93.05 

C6786 

B.  oklahomensis 

4313.1 

43.113.1 

43.113.1 

40.713.1 

100 

91.17 

91.33 

K96243 

B.  pseudomallei 

48.6  ±3.1 

48.613.2 

48.613.2 

45.6  ±3.3 

39.613.1 

100 

99.04 

ATCC23344 

B.  mallei 

48.7  ±3.1 

49.113.1 

49.113.1 

45.7  ±3.2 

40.313.1 

92.512.8 

100 

375 

376 

377 
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MALDI-TOF: 


Preparations  of  bacterial  isolates  (MSMB43,  MSMB121,  MSMB122)  for  matrix-assisted 
laser-desorption/ionization  time-of- flight  mass  spectrometry  (MALDI-TOF  MS)  were  done 
according  to  the  ethanol/fonnic  acid  extraction  protocol  recommended  by  the  manufacturer 
(Bruker  Daltonics)  and  as  previously  described  (Mellmann  et  ah,  2008)  Briefly,  a  loopful  of 
bacterial  material  was  evenly  dissolved  in  300  pi  analytical  grade  water,  and  900  pi  pure  ethanol 
was  added.  The  cell  suspension  was  centrifuged  at  13,000  x  g  for  2  minutes,  and  the  supernatant 
was  discarded.  The  centrifugation  was  repeated,  and  the  residual  ethanol  was  discarded.  The 
pellet  was  air  dried  and  thoroughly  resuspended  in  5-50  pi  70%  formic  acid  depending  on  the 
size,  and,  finally,  an  equal  volume  of  acetonitrile  was  added.  After  centrifugation  at  13,000  x  g 
for  2  minutes,  1  pi  of  the  supernatant  was  transferred  to  the  MALDI  target  plate.  After  air-drying 
at  room  temperature  1  pi  of  matrix  solution  (saturated  solution  of  a-cyano-hydroxy-cinnaminic 
acid  in  50%  aqueous  acetonitrile  containing  2.5%  trifluoroacetic  acid)  per  spot  were  was  applied. 
MALDI-TOF  MS  was  conducted  using  a  Microflex  LT  mass  spectrometer  (Bruker  Daltonics) 
equipped  with  a  nitrogen  laser.  All  spectra  were  recorded  in  linear,  positive  ion  mode  across  a 
mass/charge  ratio  (m/z)  of  2,000  to  20,000.  The  acceleration  voltage  was  20  kV.  Spectra  were 
collected  as  a  sum  of  240  shots  across  a  spot.  Main  spectra  were  calculated  from  8  spectra  per 
strain  and  used  for  construction  of  a  score  oriented  dendrogram  by  using  the  BioTyper  software 
(version  2.3,  Bruker  Daltonics).  In  these  analyses  B.  humptydooensis  sp.  nov.  strains  grouped 
together  with  B.  ubonensis  from  the  Bee  and  other  members  of  the  B.  pseudomallei  complex 
(Fig.  SI). 


Fatty  acid  methyl  ester  (FAME)  profile  analysis: 


Fatty  acid  methyl  esters  analysis  was  unable  to  distinguish  among  the  fatty  acid  profiles  from  the 
three  B.  humptydooensis  sp.  nov.  strains  and  its  closely  related  species  (five  B.  thailandensis,  two 
B.  oklahomensis,  and  three  B.  ubonensis  strains).  The  fatty  acid  methyl  esters  were  performed 
using  the  MIDI  Sherlock"  Microbial  Identification  System  by  Microbial  ID,  Inc.  (Delaware, 
USA).  Similarity  index  analysis  was  used  to  display  the  relatedness  of  FAME  profiles  produced 
by  these  species.  The  analysis  revealed  that  FAME  profiles  of  the  four  B.  humptydooensis  sp. 
nov.  strains  grouped  into  two  sub-groups,  along  with  strains  from  other  species  (Fig.  S2).  This 
finding  suggests  that  B.  humptydooensis  sp.  nov.  had  similar  FAME  profiles  with  other  tested 
genetically  related  Burkholderia  species  and  this  technique  cannot  be  used  to  distinguish  B. 
humptydooensis  sp.  nov.  or  distinguish  the  other  closely  related  species. 

In  conclusion,  we  have  utilized  comprehensive  genotyping  techniques  including  16S 
rRNA,  recA ,  MLST,  and  whole  genome  based  GGDC  to  further  support  the  existence  of  a  new 
species  that  is  distinct  but  genetically  related  to  the  four  members  of  the  pseudomallei  complex 
(B.  pseudomallei ,  B.  mallei,  B.  thailandensis,  and  B.  oklahomensis).  These  analyses  confirm  the 
speciation  of  B.  humptydooensis  sp.  nov.,  a  soil  bacterial  saprophyte  found  in  the  Northern 
Territory  of  Australia  where  melioidosis  is  highly  endemic.  The  addition  of  B.  humptydooensis 
sp.  nov.  as  a  new  member  of  the  pseudomallei  complex  will  benefit  evolutionary  studies  of  B. 
pseudomallei,  the  serious  bacterial  pathogen  that  shares  a  similar  ecological  niche  with  this  new 


species. 


Fig.  SI.  Dendrogram  demonstrating  strain  relatedness  revealed  by  MALDI  -  TOF  analysis. 
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Fig.  S2.  A  dendrogram  demonstrating  the  relatedness  of  fatty  acid  compositions  in  B. 
humptydooensis  sp.  nov.  and  other  closely  related  species.  Similarity  Index  in  the  Microbial 
Identification  System  (MIS)  was  used  to  display  how  closely  the  fatty  acid  compositions  found 
in  B.  humptydooensis  sp.  nov.,  a  new  species,  compared  with  those  from  other  genetically  related 
species.  FAME  (Fatty  Acid  Methyl  Esters)  profiles  were  generated  using  GC-MS.  The 
dendrogram  was  generated  based  upon  the  clustering  analysis  technique  to  produce  unweighted 
pair  matching  based  on  FAME  profiles.  Our  analysis  has  shown  that  FAME  profiles  from  three 
tested  B.  humptydooensis  sp.  nov.  strains  were  grouped  into  two  sub-groups.  The  first  group  also 
contained  one  B.  ubonensis  and  two  B.  oklahomensis  strains,  while  the  second  group  contained 
most  of  B.  thailandensis  strains  and  two  B.  ubonensis  strains.  We  noted  that  B.  thailandensis 
E264  had  a  distinct  FAME  profile  to  both  subgroups. 
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Fig.  S3.  16S  maximum  likelihood  (500  bootstrap)  phylogeny  using  both  copies  of  16S  from  each 
B.  humptydooensis  sp.  nov.  strain  (1420  bp  and  14  sequences). 
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Fig.  S4.  Maximum  likelihood  phylogeny  (1,500  bootstrap)  of  recA  sequence  (335  bp)  using  193 
sequences. 


